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Abstract 

We study the behaviour of the solutions of the stationary diffusion equa- 
tion as a function of a possibly rough (L y -) diffusivity. This includes the 
boundary behaviour of the solution maps, associating to each diffusivity the 
solution corresponding to some fixed source function, when the diffusivity 
approaches infinite values in parts of the medium. In n-dimensions, n ^ 1, 
by assuming a weak notion of convergence on the set of diffusivities, we 
prove the strong sequential continuity of the solution maps. In 1-dimension, 
we prove a stronger result, i.e., the unique extendability of the map of so- 
lution operators, associating to each diffusivity the corresponding solution 
operator, to a sequentially continuous map in the operator norm on a set con- 
taining 'diffusivities' assuming infinite values in parts of the medium. In 
this case, we also give explicit estimates on the convergence behaviour of 
the map. 
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1 Motivation 



Numerical methods for the diffusion equation with rough coefficients have been 
studied extensively [3, 4, 5, 12, 13, 14, 16] in the preconditioning (multigrid, 
domain decomposition, and related iterative methods) literature starting the early 
eighties and still continue to be an active area of research in various precondition- 
ing efforts [20, 21]. This article came about out of a need of deeper understanding 
of the performance of preconditioners and their connection to the underlying PDE. 

In a recent article [2], the first author constructed a new preconditioning strategy 
with rigorous justification which is comparable to algebraic multigrid. It is shown 
in [2] that analytical tools such as singular perturbation analysis gives valuable in- 
sight about the asymptotic behavior of the solution of the underlying PDE, hence, 
provides feedback for preconditioner construction. 

According to experience, the performance of a preconditioner depends essentially 
on the degree to which the preconditioner operator approximates the underlying 
operator. Then, the fundamental need is to explain the effectiveness of the pre- 
conditioner and to justify that rigorously. In that respect, one can view the tools 
in this article as steps towards adding tools to the arsenal of methods of analysis 
for rigorous justification at the interface of preconditioning and operator theory. 
Direct connections from the results here to preconditioning will be the subject for 
future research. 



2 Introduction 

The diffusion equation 

— =div(pgradu)+/ (2.0.1) 

describes general diffusion processes, including the propagation of heat, and flows 
through porous media. Here u is the density of the diffusing material, p is the 
diffusivity of the material, and the function / describes the distribution of 'sources' 
and 'sinks'. This paper focuses on stationary solutions of (2.0.1) satisfying 

— div (pgrad u) = f . (2.0.2) 

For instance, the fictitious domain method and composite materials are sources of 
rough coefficients; see the references in [14]. Important current applications deal 
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with composite materials whose components have nearly constant diffusivity, but 
vary by several orders of magnitude. In composite material applications, it is quite 
common to idealize the diffusivity by a piecewise constant function and also to 
consider limits where the values of that function approach zero or infinity in parts 
of the material. 

Results of such study were given first by J. L. Lions [15]. In his lecture notes, 
he considers the limit of the solution of (2.0.2) where the limit is associated to 
a one-parameter family of piecewise constant diffusivities approaching zero on 
a subdomain. The same piecewise constant one-parametric approach was used 
in [4, 13], but with diffusivities approaching infinity on a subdomain. The limi- 
tation of the one-parametric approach is its dependence on the particular approxi- 
mating sequence. To the knowledge of the authors, this paper is the first to address 
these questions in the necessary generality. Hence, we consider general fami- 
lies of diffusivities that are not necessarily piecewise constant. In addition, due 
to the atomistic structure of matter, the physical treatment of diffusion involves 
regular (C 00 -) diffusivity. It is unclear to what extent the idealization of diffusiv- 
ity by piecewise constant coefficients has the capability to capture the underlying 
physics. Mathematically, the severe contrast in diffusivity should be represented 
by a regular function whose size is changing drastically over small distances in 
interface regions. In this paper, we demonstrate that the assumption of piecewise 
constant diffusivities is meaningful by showing a continuous dependence of the 
solutions on the diffusivity. 

Furthermore, the diffusion equation is meaningless if the 'diffusivity' is infinite 
or zero in regions of the material. Physics requires nowhere vanishing diffusivity 
in the interior of the material. As a consequence, only the relative size of diffusiv- 
ities should be significant. Therefore, physically, one might expect that both types 
of the above limits are equivalent, but mathematically there are differences. The 
limit of the solution as diffusivity approaches infinite values exists. However, only 
the limit of the scaled solution exists as diffusivity approaches zero values (see 
Example 4.0.1 1 for both cases). That is why, we choose to work with diffusivity 
approaching infinity. We will refer these cases as 'asymptotic cases'. 

Also, the treatment in [4, 13] considers only limits on specific parts of the ma- 
terial. In this connection, it should also be remarked that, although (2.0.1), (2.0.2) 
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are linear equations, in general, their solutions depend non-linearly on the coeffi- 
cients. 

For the treatment of these questions, we use methods from operator theory. For 
this, we use a common approach to give (2.0.1) a well-defined meaning that, in a 
first step, represents the diffusion operator 

— divpgrad (2.0.3) 

as a densely-defined positive self-adjoint linear operator A p in a suitable Hilbert 
space. As a result, (2.0.2) is represented by the equation 

A p u = f , 

where / is an element of the Hilbert space, and u is from the domain, D(A P ), of 
A 1 

Specifically, we treat the class C of diffusivities p e L°°(f2) that are almost ev- 
erywhere ^ eonfi for some e > 0, where Q a M. n , n e N*, is some non-empty 
open subset. By use of Dirichlet boundary conditions, it defines A p as an oper- 
ator in the complex Hilbert space L^(Q). For non-smooth p, the domain of A p 
depends heavily on p. This fact significantly complicates the study of sequences 
of functions of A p . 

In this paper, we turn to a first-order formulation of (2.0.2) which is often referred 
as mixed formulation in the discretization literature [7]. The first-order formula- 
tion was popularized in the least squares finite element community by the so-called 
FOSLS pioneering paper [9]. Here, we provide the self-adjointness of the corre- 
sponding operator A p in a Hilbert space. The key property of A p is that its domain, 
D(Ap), is independent of p. This property is exploited in establishing the conti- 
nuity of the solutions A~ l f as a function of p. Moreover, A p remains defined for 
the asymptotic cases when (2.0.1), (2.0.2) are ill-defined. This fact is used in the 
study of the asymptotic cases. 



After that, the abstract theory of strongly continuous one-parameter semigroups of operators can 
be used to associate a rigorous formulation of a well-posed initial value problem to (2.0.1) [6, 10, 
17]. In this, A p becomes the infinitesimal generator of time evolution. This last step will not be 
detailed here. 
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Specifically, for p e C and by assuming a weak notion of convergence in C, 
we show that the maps that associate p to the operator A~ l and — pVA^ 1 , re- 
spectively, are strongly sequentially continuous, see Theorem 5.0.21 and Corol- 
lary 5.0.22. In particular, this shows in these cases that the approximation by dis- 
continuous coefficients to physical diffusivity is indeed meaningful. In addition, 
for the case n = 1 and bounded open intervals of M, we show stronger results that 
include also the asymptotic cases, except that where the asymptotic 'diffusivity' 
is almost everywhere infinite on the interval. In this case, the maps that associate 
p to the operator Ay- and — (1/p) VAy\, respectively, have unique extensions 
to sequentially continuous maps in the operator norm on the set of a.e. positive 
elements of L°°(f2)\{0}, see Corollary 6.0.26, 6.0.27. In addition, an explicit es- 
timate of the convergence behaviour of the maps is given, see Theorem 6.0.25. It 
is still an open problem, whether the last results are generalizable to dimensions 
n 5s 2. 

3 Basic notation 

Mainly, this section introduces basic notation. In particular, an operator theoretic 
definition of Sobolev spaces is given that is based on weak derivative operators, 
instead of distributions. In such formulation, the completeness of the Sobolev 
spaces is an obvious consequence of the cussedness of these operators. Also, we 
give some basic results that are connected to this formulation. For the convenience 
to the reader, corresponding proofs are given in the appendix. 

General Assumption 3.0.1. In the following, let n e N* and Q be a non-empty 
open subset of R n . 

We follow common usage and do not differentiate between a function / which is 
almost everywhere defined (with respect to a chosen measure) on some set and 
the associated equivalence class consisting of all functions which are almost ev- 
erywhere defined on that set and differ from / only on a set of measure zero. The 
following definitions need to be understood in this sense. 

Definition 3.0.2. (Complex L p -spaces) 

(i) For p > 0, the symbol L^(Q) denotes the vector space of all complex- 
valued measurable functions / which are a.e. defined on SI and such that 
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\f\ p is integrable with respect to the Lebesgue measure v n . For every such 
/, we define the L p -norm |/| p corresponding to / by 



up •" 



f l/l 
Jn 



i/p 

p dv n " 



In addition, for the special case p = 2, we define a scalar product ( | ) 2 on 
4(fi)by 

(f\g) 2 -= f f*gdv n , 



for all f,g e L^(fi). Here * denotes complex conjugation on C. As a 
consequence, ( | ) 2 is antilinear in the first argument and linear in its second. 
This convention will be used for sesquilinear forms in general. 

(ii) Lj?(f2) denotes the vector space of complex-valued measurable bounded 
functions on Q. For every / e Lq(Q), we define 

|/|co:= sup 1/(3)1 . 
xefl 

(iii) For every k e N* and f,ge (L^(f2)) fe , we define 

k / k \ 

(f\9Xk--=t<fi\9i)2 ■ Uhk--= Sll/ill! • 
i=i \j'=i / 

Definition 3.0.3. (Weak derivatives and Sobolev spaces) We define 

(i) for every multi-index a e N n the densely-defined linear operator 5 a in 

4(n)by 

where * denotes the adjoint operation and 

n 

\a\ := J] a j ■ 
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(iii) for every k e N the Sobolev space W^(Q) of order k by 

W^) := P| D(d a ) . 

aeN n ,\a\^k 

Equipped with the scalar product 

<,> fc : (W^)f^C, 

defined by 

</,4 := J] (d a f\d a g) 2 

aeN",|«|s; k 

for all f,g 6 W^(J1), W^(n) becomes a Hilbert space. 

(iv) Wq C (Q) as the closure of C™(0, C) in (W^(Q), \\\ \\\k), where ||| ||| fc denotes 
the norm that is induced on W£(Q) by ( , ) k . 

We note that 

Lemma 3.0.4. (Partial integration) 

{f\d ek g\ = -{d ek f\g\ (3.0.4) 

for all (/, g) e Wq c (fi) x Wq(Q) and k 6 N*, where denotes the /c-th canonical 
unit vector of M n . 

The next defines gradient operators. 

Definition 3.0.5. (Gradient operators) We define the (L^(r2)) n -valued densely- 
defined linear operators in £^(0) 

V : C °°(ft,C) - (L 2 c (n)) n , V w : W&Sl) -» (4(0))" 

by 

vo/ := (H? ■ • • ' J£) ' Vw5 := t(aei5 ' ■ ■ ■ ' 

for all / e qf (O, C) and g e W^(Q). 

Then the following holds. 

Lemma 3.0.6. (Adjoints of gradient operators) 

(V *)* = Vj ! rm , Vj i ,o) = V * . (3.0.5) 
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4 Basic properties of the diffusion operator 

This section provides the basis of the paper. It defines the diffusion operator as 
operator in Lq(Q) and gives basic properties. 

Definition 4.0.7. Let p : 17 — > R be measurable and such that 1 jp is a.e. defined 
on fl We define the linear operator A : D(A) — > Lq(Q) in Lg(Q) by 

D(A):= {u e WjcCfl) = (l/p)V„u eC(V *)} 

and 

Au:= V *(l/p) V„,u 

for every u 6 -D(A). 

Diffusion operators corresponding to diffusivities from the following large subset 
£ of L°°(fi) will turn out to be densely-defined self-adjoint linear operators. 

Definition 4.0.8. We define the subset C of L°°(fi) to consist of those elements p 
for which there are real C\ , C2 satisfying C2 ^ C\ > and such that C\ ^ p ^ 
C2 a.e. on fl Note that the last also implies that 1/p e C and in particular that 

1/C 2 < l/p< 1/Ci a.e. on 0. 

The next proves the self-adjointness of diffusion operators corresponding to diffu- 
sivities from C. For this, so called 'form methods' from operator theory are used. 
For these methods, see [11]. 

Theorem 4.0.9. Let p e C. Then A is a densely-defined linear self-adjoint opera- 
tor in Ljc(n). 

Proof. For this, we define a positive Hermitian sesquilinear form s : (Wq C (CI)) 2 — 
Cby 

s(u, v) := (V w u I (1/p) V^t;> 2in 
for all u,vs Wl c (n). Then < | ) s : (Wl c (n)) 2 -> C, defined by 

<u|w> s := +<n|u> 2 

for every u, v s Wg c (f2), defines a scalar product on Wq C (Q) with induced norm 

|| || s : Wq )C (0) — > R given by 

\\ u fs = i^wu I (1/p) y w u\ n + \\u\\l 
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for all u 6 W c (fi). In particular, s is closable. For the proof, let Ui, U2, ■ ■ ■ be a 
Cauchy sequence in (Wq c (f2), || | s ) and such that 

lim \\u u \\2 = . 

v— >oo 

We note that 

min{l, 1/C 2 } •< 77- ||V^n|| 2 ,n + H2 < Mis 

< 7T ||V w «||2,n + Mil ^ max{l, 1/Ci} fuf \ , 

where Ci, C2 e! satisfy C2 ^ Ci > and are such that Ci ^ p ^ C2 a.e. on Q,, 
and hence that j| | s and the restriction of ||| to W 1 c (O) are equivalent. Hence 
it follows that 

lim lluJL = . 

U—KX) 

Since (Wq c (fi), || || s ) is in particular complete, it follows that s coincides with its 
closure. As a consequence, there is a unique densely-defined linear self-adjoint 
operator A : D(A) —> L^(Q) in L^(Q) such that D (A) is a dense subspace of 
(Wq C {Q), I |||i) and such that 

(u\Au) 2 = s(u,u) = (V w u I (1/p) V w u) 2 n 

for all u e D(A). In particular, D(A) consists of all u e Wq c (f2) for which there 
is / e L^CD.) such that 

(/ I • • • >2 llVg C (f2) = ( (V^) V ™ M I Vu, . . . >2,n|vy ( i c (f7) ' 

Further, if u and / satisfy these requirements, then 

Au = f . 

Hence u e D(A) if and only if 

(1/p) V w « e D ((v.|^ c(n) ) ) = D(V *) 

and in this case 

Au = V *(l/p)V w u . 

□ 
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Fig. 1: Graphs of u from Example 4.0.1 1 as a function of M, 



For completeness, the next gives the proof that diffusion operators corresponding 
to diffusivities from C have a purely discrete spectrum, i.e., that their spectrum is a 
discrete subset of the real numbers consisting of eigenvalues of finite multiplicity 
and that there is a Hilbert basis consisting of eigenvectors. This result is not used 
in the following. 

Corollary 4.0.10. Let p e £ and, in addition, 17 be bounded. Then A has a purely 
discrete spectrum. 

Proof. According to the proof of Theorem 4.0.9, || || s : Wq c (f2) — > R defines a 
norm which is equivalent to the restriction of ||| |j|i to Wq C (U). Hence the closed 
unit ball B in (Wq C (£1), \\ \\ s ) is contained in a closed ball of (Wq C (£1), ||| |||i). 
The last is relatively compact in L^(J7). From this, it follows also that B is rel- 
atively compact in L^(il). Hence it follows, see, e.g., [18] Vol. IV, that A has a 
purely discrete spectrum. □ 

Example 4.0.11. The following example illustrates the influence of discontinuities 
of the diffusivity on the regularity of the elements in D(A). Consider the case that 
= /:=(— 1,1) and a piecewise constant diffusivity p : I — » M given by 



1 



if -1 < x < -1/2 
if -1/2 ^ x ^ 1/2 
if 1/2 < x < 1 



p(x) : 



= < 



M 



1 
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for x £ I, where M > 0. Then Au = f, where u : I — > R is defined by 

f(l-x 2 )/2 if -1 < x s$ -1/2 

u(x) : = J (l _ Ax 2 + 3M)/(8M) if -1/2 < x < 1/2 
[(l-x 2 )/2 ifl/2<o;<l 

and / is the constant function on I of value 1. We note that u' has no extension to 
a continuous function on i" if M ^ 1. In general, discontinuities in the diffusivity 
cause low regularity of elements in D(A). Also, see the concluding remarks. 

There is a unique solution Uf to the equation 

Au f = f 

for every / e L%*(Q) if and only if A is bijective or equivalently, if and only if 
is not part of the spectrum of A. In general, A is not bijective. For instance, the 
operator A that is associated to = W 1 and the diffusivity p(x) = 1 for every 
x e M n is not surjective. Below, we place a restriction on Q, that leads to bijective 
diffusion operators. 

General Assumption 4.0.12. In the following, we assume that Q is in addition 
such that the following Poincare inequality is valid 

\\d^f\\ 2 ^c\\f\\ 2 (4.0.6) 

for some j £ {1, . . . , n} and every / £ Wq C (Q), where c > 0. In the remainder, 
such c is considered chosen. 

Remark 4.0.13. It is known that SI of the assumed type are not necessarily bounded. 
For instance, every non-trivial open set, for which there is n £ M™\{0} along with 
real numbers a, b such that 

a < x ■ n < b 

for all x £ 0, is of this type. 

In particular, the following proves that diffusion operators corresponding to diffu- 
sivities from C are bijective. 

Theorem 4.0.14. Let p £ C. The spectrum a (A) of A satisfies 

a {A) c [c 2 /C,oo) , (4.0.7) 
where j £ {1, . . . , n} is such that rij ^ and C > is such that p ^ C a.e. on Q. 
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Proof. For this, let j e {1, . . . , n} be such that rij ^ 0. For u e D(A), it follows 
that 

(u\Au) 2 = (V w u | (1/p) W w u) 2n ^ C- 1 \\V w u ||| n ^ c 2 C" 1 ||u|| , 

where C > is such that p ^ C a.e. on 17. Hence it follows the validity of 
(4.0.7). □ 

5 Properties of a first order operator connected to the dif- 
fusion operator 

As indicated by Example 4.0. 1 1 , for non-smooth diffusivities p, the condition that 
pV w u e Wji(O) in the definition of the domain of A leads to a strong dependence 
of that domain on the diffusivity. This fact poses an obstacle to the study of the 
map, associating to every diffusivity p e C the corresponding operator A^ 1 , by 
the notion of strong resolvent convergence, see, [18, Volume I, Section VIII.7], 
[11, Section VIII, §1]. By use of the following vector partial differential operator 
of the first order A, this problem can be circumvented. Its domain is indepen- 
dent of the diffusivity. The connection of the resolvents of A and A is given in 
Theorem 5.0.20. 

Definition 5.0.15. Let p e L 00 ^). We define the densely-defined, linear operator 

A : WjcCfi) x D(V *) - Ll(Q) x (L 2 : (f2)) n in L 2 («) x (4(0))" by 

A(u,q) : = ( V *<?, V w u-pq) 
for every (u,q) e Wo,c( n ) x ^( v o*)- 
Theorem 5.0.16. The operator A is self-adjoint. 

Proof. The statement is a consequence of Lemma 3.0.6. □ 

The following gives a characterization of the kernel of A. In particular, the result 
implies that A is bijective for diffusivities from C. 

Theorem 5.0.17. Letp e L c0 (il) be a.e. positive. Then 

keri = {0} x (ker V * n kerT? ) , 

where Tp e L(L 2 (r2), L^(fi)) denotes the maximal multiplication operator in 
tnat i s associated to p. 
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Proof. *c': Let q 6 ker V * n ker Then (0, g) e £>(A) and 

V *g = 0, ~pg = 0. 

Hence it follows that (0, g) e ker A. 
'=>': Let (u, q) £ ker A Then 

V *q = 0,V w u-pq = (5.0.8) 

and hence 

= <g| V w u-pq\ n = (q\V w u\ n -(q\pq) 2in 

= < V *g | M>2 - ||p 1/2 <? || 2 ,n = -||p 1/2 g ||2,n • 

The last implies that 

gekerTl 1 ^ , 

where T^i/a 6 L(L 2 : (J7), L^(fi)) denotes the maximal multiplication operator in 
Ifc(^) that is associated to p 1 / 2 , and hence also that 

gekerT? . 

Further, by (5.0.8)-2), it follows that 

V w u = . 

The last implies that 

V *V„,u = 

and hence by Theorem 4.0. 14 that u = 0. □ 

The following example shows that the kernel of A is non-trivial if p vanishes on 
some open subset of Q,. The vanishing of p on non-empty subsets of Q, corresponds 
to the asymptotic cases mentioned in the introduction. 

Example 5.0.18. In the following, we give q e C^°(R n ,R") n ker V * for n 5s 2. 
For this, let h be an element of Cg°(M) with support contained in [—1, 1]. In addi- 
tion, let a be a non-zero antisymmetric n x n-matrix. We define q e C^°(IR n , M. n ) 
by 

, , MM 2 ) v 



13 



for all x = (xi, . . . , x n ) 6 W 1 . Then 

i=l % i,j = l 

for all x = (xx, . . . , x n ) e W 1 and hence q e ker Vo*. 

The following lemma prepares the subsequent theorem which estimates the size of 
the gap around in the spectrum of A and gives a representation of the resolvent 
of A in terms of the resolvent of A, i.e., (5.0.9). The main tool in the proof is the 
Closed Graph Theorem in the form of see Theorem 3.1.9 in [6]. 

Lemma 5.0.19. Let e C, u(A) the spectrum of A, A < mm{a(A)} and A\ : = 
A - A. Then 

(i) V w A?eL(Ll(n),(L*(n)r), 

(ii) A^V * = (V^ 1 )*, 

(iii) D(A\ /2 ) = WjcCfi) and A 1 / 2 : W£ C (Q) -> L^(Q) is continuous, 

(iv) VwA^Vo* is a positive self-adjoint element of L((L 2 c (n)) n , (L£(J2)) n ). 

Proof. '({)': Since A e C\a(A), A\ is densely-defined, linear and bijective. 
Further, Ax is self-adjoint and strictly positive. As a consequence of its self- 
adjointness, A\ is in particular closed. Further, according to Lemma 3.0.6, the 
restriction of V w to Wq c (fi) is a closed linear operator in L^(Q) with values in 
{L^(p.)) n . Since D(A) a Wg C (fi) it follows from the closed graph theorem, 
e.g., see Theorem 3.16 in [6], the existence of C e [0, oo) such that 

\\V w f\\ 2 ,n^C\\A X f\\ 2 

for all / e D(A). As a consequence, it follows for every / 6 Lq(£1) that 

\\V w A- 1 f\\ 2<n ^C\\A x A- 1 f\\ 2 = C\\f\\ 2 
and hence that V^A" 1 e L{I? c (ti), (i^(fi)) n ). 

'(ii)': A^Vq* is a densely-defined, linear operator in (L^(f2))™ with values in 
Ll(n). Further, it follows for / 6 L&(fi) that 

< / | A- 1 V *g > 2 = < A- 1 / | V *<? > 2 = < V^ A V | q \ n 
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for every q e £>(V *) and hence that f e D^A^Vq*)*) as well as that 
As a consequence, 

(A^Vo*)* = . 
In particular, AT 1 Vq* is closable and 



A^v * = {v w A- l y . 

'(iii)': In a first step, we prove the statement for the case A = 0. For this, we note 
that, as a consequence of Theorem 4.0.14, < mm{cr(74)}. Further, we note that 
D(A) is a core A 1 / 2 . For instance, this follows by Theorem 3.1.9 in [6]. Hence 
D(A) is dense in the Banach space (D(A 1 ^ 2 ), || H^i/2), where 

ll/IUv 2 := [11/111 + \\A l / 2 ff 2 f /2 

for every / 6 D(A 1 / 2 ). Further, it follows for / e D(A) that 

U 1/2 f\\l = (f\Af) 2 = (V w f\ (1/p) V w f \ n = s(f, f) , 

where the real numbers C\ , C 2 and the sesquilinear form s are as in the proof of 
Theorem 4.0.9, and hence that 

min{l,l/C 2 } HI/HI? < \\ff Al/2 ^ma^MM} Ill/Ill? • 

As a consequence, the restrictions of || || A i/2 and ||| to D(A) are equivalent. 
Since D(A) is dense in (D(A 1 I 2 ), || || A i /2 ), it follows for / e ^(A 1 / 2 ) the exis- 
tence of a sequence fi, f 2 , ■ ■ ■ in D(A) such that 

Jim ||/,-/|| A i/2 =0. 

Since the inclusion of (D(A 1 / 2 ), || H^i/2) into -L^(^) * s continuous, this implies 
also that 

Jim ||/„~/! 2 = 0. 

Since the restrictions of || ||^i/2 and ||| to D(A) are equivalent, it follows that 
fx, f 2 , ■ ■ ■ is a Cauchy sequence in Wq C (Q) and hence convergent to some / e 
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Wq(0). Since the embedding of (W^(fi), ||| |||i) into L^(£l) is continuous, it 
follows also that 

lim WU-fh =0 

>co 

and hence that / = / 6 Wo(f2). Further, it follows that 

min{l,l/C 2 } HI/I? ^ \\ff A y 2 ^max{l,l/Ci} |||/|||? 

and hence that || || A i/2 and the restriction of ||| \ 1 to D(A 1 / 2 ) are equivalent. Since 
according to the proof of Theorem 4.0.9, D(A) is a dense subspace of (W 1 C (Q), 
I |||i), we conclude that Z^A 1 / 2 ) = Wjcffi) and that A 1 / 2 : W^ C (Q) -> L^(ft) 
is continuous. From this, we conclude that statement of (ii) as follows. For this, let 
A e R\<t(A) such that A > max{0, A}. Since M\cr(.A) is open, such A exists. We 
note that D(A) is a core also for A 1 ^ 2 and A 1 / 2 - For instance, this follows by The- 



orem 3.1.9 in [6]. Hence D(A) is dense in the Banach spaces (D(A 
P(^l /2 ),|| L1/2), where 



1/2 X 

A )i\\\\ A V 2 )> 



.1/2 



[]|/ll + ]l4 /2 /lli] 1/2 , M A * ^[Ml + WAfgll] 1 ' 2 , 

A 



for all / 6 D(A\ /2 ) and g e D(a][ 2 ). Further, it follows for every / e D(A) that 

11/11^ = ii4 /2 /iii + mi = (f\A X f\ + n/ii 

= </|A A /> 2 + ||/| + (A - A)||/||| = ||/||^ /2 + (A - A)||/| 
and hence that 



as well as that 



r 1/2 ^ [i + (a - a)] h/ii^ 1/2 



Since D(A) is dense in the Banach spaces (D(A 1 ^ 2 ), \\ || 1/2) and (D(A 1 / 2 ), \\ \\ 1/2] 

A A 

it follows that 

D(Af) = DiA 1 / 2 ) 
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D(A\ /2 ) = D{A l l 2 ) 



as well as the equivalence of the norms II II i/ 2 and II II .1/2. In particular, this im- 
plies that 

and the equivalence of the norms || || 1/2 and || l^i/a. By this, the statement (ii) 

A 

follows from the corresponding statement of (ii) for the special case that A = 0. 
'(iv)': In a first step, we conclude that 

V W A- 1/2 e L(4(tt),(4(0)n . 

As a consequence of the analogous properties of A\, A 1 ^ 2 is densely-defined, 

linear, self-adjoint and bijective. Since A 1 ^ 2 is self-adjoint, A 1 ^ 2 is in particular- 
closed. Further, according to Lemma 3.0.6, the restriction of V w to WQ C (fi) is 

a closed linear operator in L^(U) with values in (L^O))™. Since D{A\[ 2 ) = 
Wq C (Q), it follows from the closed graph theorem, e.g., see Theorem 3.16 in [6], 
the existence of C e [0, 00) such that 

||V w /|| 2 , n ^C||4 /2 /|| 2 

for all / 6 D(A 1 ^ 2 ). As a consequence, it follows for every / e L^(fi) that 

||V^ 1/2 /|| 2 ,n ^CWAfA-^fW, = C\\f\\ 2 

and hence that V^A^ 1 ^ 2 e L(L^(fl), (X^fi))"). In a second step, we conclude 
that 

A- x 1/2 V *eL((L 2 c (n)r,L 2 c (n)) . 

A x 1//2 Vo* is a densely-defined, linear operator in (L|.(0)) n with values in L^(Q). 
Further, it follows for / e L^(Q) that 

< / I A~ 1/2 V *q > 2 = < A~ 1/2 f \V *q) 2 = < V w A~ 1/2 f I q \ n 
for every qeD(V Q *) and hence that / e L>((A~ 1/2 V *)*) as well as that 



As a consequence, 



(^- 1/2 V *)*/ = V w A- 1/2 f 



(A" 1/2 V *)* = V W A~ X 1/2 
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In particular, A x 1 ^ 2 Vq* is closable and 



A-VVo* = ( V^" 1/2 )* e L((Ll(n)r, L&Sl)) . 
Further, we note that 

V^Vo*/ = V^ 1/2 ^ 1/2 V 7 = V w A- l/2 (V w A~ 1/2 )*f 



for every / e D(\7 *). Hence it follows that V W A X 1 Vo* is a positive self-adjoint 

£(n))»,(z£( 



element of L((L£ (ft)) n , (1^(0))™) . □ 



By help of the previous lemma, we can now estimate the size of the gap around 
in the spectrum of A and give a representation of the resolvent of A in terms of the 
resolvent of A, i.e., (5.0.9). 

Theorem 5.0.20. Let p e C, Ci, C 2 6 R satisfy C 2 ^ Ci > and be such that 
C\ ^ p < C 2 a.e. on 0. Further, let j e {1, . . . , n} be such that rij ^ 0. Then the 
interval 

J:=(-C 1 ,c 2 /(c + C 2 )) 

is contained in the resolvent set of A. In particular for X e J, (A — A) -1 is given 
by 

(A-\)-\f,g) = 

((A px - A)- 1 / + (A Px - Xy^pxg, (5.0.9) 
-Px9 +P\V w (A Px - A)" 1 / +p x y w (A Px - A)- 1 V *p A 5) 

for all (f,g) e -L^(O) x (L^(fi)) n , where p x ■= p + \, Px ■= l/p x , and A Px is 
the operator corresponding to p\ according to Definition 4.0.7. 

Proof. For this, let A 6 J. Then, 

0<A + Ci^p + A^A + C 2 

and p\ := p + A 6 C. Further, we denote by A Px the operator corresponding to p\ 
according to Definition 4.0.7. As a consequence of Theorem 4.0.14, the spectrum 
of A Px — A is contained in the interval 

[c 2 (A + C 2 )- 1 - A,oo) , 
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The inequality 

c 2 (\ + c 2 y 1 - A> 

is equivalent to 

(A + f) 2 -f = A(A + C 2 )<c 2 . 
The last is equivalent to 

c2 + f-T <A< V c + f"T- 

We note that 
and that 

Hence it follows that A Px — A is self-adjoint, strictly positive and bijective. We 
define the bounded linear operator B e L(L 2 c (n) x (Ll(Q)) n , L 2 c (n) x (L 2 z (p)) n ) 
by 

B(f,g) := {{A px - A)" 1 / + (A Px - p\g , 

- P\9 + w {A Px - \y l f+P\V w {A Px - A)- 1 V *pa^) 

for all (/,c/) e L^(ft) x (L^(0)) n , where p A := l/p A 6 L^O). Further, we 
define the subspace D of L^ft) x (L^(0)) n by 

L> :={(/, 9) 6 4(^) x = 6 ^(V *)} . 

We note that the subspace 

{p x g:geC?(n,C)} 

of L^(Q) is dense in L^(Q). For the proof, let / 6 Lf^(U). Since p x e L G0 (fi), 
j»a/ ^ ^c(^)- Further, since C^(Q, C) is dense in L^(fi), there exists a sequence 
/1, / 2 , . . . in C) such that 

lim \\f v -p x fh = • 
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Since for every i/eN* 

\\PxU - fh = \\Px(fu-Pxf)h ^ \\PxU\\U -Pxfh , 
it follows that 

lim \\pxU - f\\ 2 = . 

v— >oo 

Hence it follows also that D is dense in L^(Q) x (L^(S7)) n . Further, for (/, g) e 
D, it follows that B(f, g) e D{A) and that 

V * [ -Pxg + PxV w {A Px -\)- 1 f+ PxV w (A Px - A)-iVo*p A ff ] 
- X[(A Px - A)- 1 / + (A Px - ApW pxg] 
= -V *pa5 + A Px (A Px - A)" 1 / + A Px (A Px - Xy^o* Px g 
- X(A Px - A)" 1 / - X(A Px - Xy^Pxg] = f 

and that 

V w [ (A Px - A)" 1 / + (A Px - A)-iV * Px g] 
-Px[-Pxg +PxV w (A Px - A)" 1 / +pxV w {A Px -Xj-^pxg] 
= V w (A Px - A)- 1 / + V„,(A P - A - A)" 1 Vo*PA5 

5 - V„,(A P - A - A)" 1 / - V^p, - A)- 1 V *pa5 = g ■ 

Hence it follows that 

(A-X)B(f,g) = (f,g) . (5.0.10) 

Further, since D is dense in L 2 C (Q) x {I? c (Vt)) n , for (/, g) e L£(fi) x (L 2 c {tt)) n , 
there is a sequence {fi, gi), (f2, g2), ■ ■ ■ in -D that is convergent to (/, g). Since £> 
is bounded, the corresponding sequence B(f\,g{), B(f 2^2), ■ ■ • is convergent to 
B(f\ g). Since A is in particular closed, it follows that B(f, g) e D(A) and that 

(A-X)B(f,g) = (f,g) . (5.0.11) 

Therefore, A — X is surjective and hence also bijective. □ 

By help of the previous theorem, the next result follows by application of a well- 
known criterion for the strong resolvent convergence of sequences of self-adjoint 
operators. 
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Theorem 5.0.21. Let pi,p2, ■ ■ ■ be a uniformly bounded sequence in C for which 
there is e > such that p u ^ e for all veM* and which converges a.e. pointwise 
on O to poo 6 C. In addition, let A± , A2 , . . . be the associated sequence of self- 
adjoint linear operators and ^oo be the self-adjoint linear operator associated to 
Poo. Then 

s- lim A' 1 = A' 1 . 

>oo 

Proof. By application of Lebesgue's dominated convergence theorem, it follows 
that 

lim \\A v (u,q) - A x (u,q) ]| = , 
u— >oo 

for all (u,q) e L^(Q) x (L^(f2))™, where || || denotes the norm on L^(Q) x 
(L^(Q)) n . From this, the statement follows from a well-known criterion for strong 
resolvent convergence of a sequence of self-adjoint linear operators, e.g., see part 
(i) of Theorem 9. 1 6 in [ 1 9] . □ 

Corollary 5.0.22. Let pi,p~2> • • ■ be a uniformly bounded sequence in C for which 
there is e > such that p v ^ e for all i/eM* and which converges a.e. pointwise 
on O to poo 6 C. In addition, let A±, A2, ■ ■ ■ be the associated sequence of self- 
adjoint linear operators and A^ be the self-adjoint linear operator associated to 
poo. Finally, let / 6 L^(n). Then 

lim {{A- 1 / - A-Vlla = lim \\ Pu V w A- l f - Poo V w A x l f\\ 2in = , (5.0.12) 

v — ►oo v— >oo 

where p u := 1/pj, for every i/eN* andpoo := l/pa> 
Proof. By Theorem 5.0.21, it follows that 

limA- 1 (f,0)=A-\f,0) , 

where A\, A2, . . . is the associated sequence of self-adjoint linear operators to 
pi , p2 , . . . and Aco is the self-adjoint linear operator associated to p. Hence (5.0. 12) 
follows by Theorem 5.0.20. □ 

6 The one-dimensional case 

In the special cases that is given by a non-empty bounded open interval of R, 
A^ 1 can be explicitly calculated. This is somewhat surprising since in this case 
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the corresponding A is a Sturm-Liouville operator and the standard method of 
calculating its inverse, e.g., see Theorem 8.26 in [19], seems not directly applicable 
for general p e C. Still, by a direct calculation of A -1 , one can give an explicit 
expression of A^ 1 by using (5.0.9). 

Theorem 6.0.23. Let a, b e R such that a < b and fi := I := (a,b). Further, 
let p e L°°(fi)\{0} be a.e. positive. Then A is bijective and has a purely discrete 
spectrum. In particular, Ar x is given by 

for every (f,g) s (L^.(I)) 2 , where 



u(x) 



rx rx |~ rx 

g(y)dy + p{ 

Ja Ja \_Jy 



u) du 

b r rb 



f(y)dy 



q(x) 



rx ( r or rb rb 

1 P(y)dy < p(x)dx f(y)dy- g{y)dy 

Ja {.Ja \_Jy J Ja 

Ja 



y 

f f(y)dy + \\P\\i 1 \f\( b p( 

J a wo L Jv 



(x) dx 



f(y)dy 



for every x s I. Also, A 1 satisfies 



\\A~ 



2(b-a) 1^1^(1 + Ml? 



Proof. For this, we define the derivative operator 

Dr:C?{I, C)^Ll(I) 

by Djf := /' for every / s C °°(/, C). In a first step, we prove an auxiliary 
result. For this, let / £ L^(I) and he C(I, C) be defined by 



h(x] 



Ja 



dy 



for every x 6 I. Further, let ip e Cq ) (I,C). By Fubini's theorem and change of 
variables, it follows that 



(h\D lV ) 2 = {h\^') 2 = ( b h*(x)<p'(x)dx= f [ I* <p'(x)f*(y)dy 

Ja Ja L Ja 



dx 
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{(x,y)eM. 2 :a^x^b ao^j/^i} 



{(x,j/)£lR 2 :asC?/sc6Aj/^xsc6} 



<p'(x)f*(y)dxdy 
ip'(x)f*(y)dxdy 



f'{x)f*{y)dx 



and hence that 



dy 



\ v{y)t{y)dy = -UW\ 

Ja 



heWc(I) and D}h = -f . 



With the help of the previous auxiliary result, we proceed in the proof of the 
lemma. Forthis, we define for every (f,g) e (L^(J)) 2 , a corresponding B(f,g) = 
(u, q) by 



q(x) := q (x) + c , u(x) := [g(y) + p(y)(q (y) + c)] 

Ja 



dy 



where 



qo(x) 



rx rb 

f{y)dy , c -.= -\\PWi 1 [g(y) +P(y)go(y)]dy 

Ja Ja 



for every x e I. By help of the auxiliary result above, it follows that (u, q) e 
(W£(J) n C(I,C)) x £)(£>;) and that 



£>f g = / , -D*u -pq = g + pq - pq = g . 



In addition, 



dy 



ub=\ [g{y) +p(y)(go(y) + c)] 

Ja 

rb rb 

(g(y) +P(y)%(y))dy + c \ p(y)dy = o. 

Ja Ja 

As a consequence, 

u a = u b = . 

From the last, it follows also that u e W 1 C (I). For the proof, let ip e C°°(M) be 
such that 



¥>((-oo,0]) c: {0} , ¥>([!, oo)) a {1} , Ran^c [0,1] . 
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Such a function is easy to construct. For i/eN*,we define ip v e C° 

<p v (x) '■= f( v{x — a)(b — x) — 1 ) 

for every x s I. Then it follows for every !/eN* satisfying v ^ b — a and x e I 
that 



if x e J \ (o + v 



-2 



.-2-1 



and hence that ip v e y± , 



1 if (x-a)(6-ac) ^ 2i/ _1 
as well as Ran(y? v ) c [0,1]. In particular, 

|(x-o)(6-a;)^(x)| 

^ 3u(a + b) (x - a)(b - x) ■ \ tp'(v(x - a)(b - x) - 1 ) | 
sj 3i/ (a + b) (x - a )(6 - s ) ■ ||y?'||oo • X {<e6/i(ie _ a)( i_« )<a/l/ , (a) 
^ 6 (a + 6) ■ ll^'lloo ■ X^f,-.^-.)^, (») 

for all x 6 /. An application of Lebesgue's dominated convergence theorem leads 
to 

lim \\(p v u— -ulU = lim \\Dj(p v u — DfulU = lim \\ip v DjU + (p'u — DfttlU = . 

^— >co >oo 1/— >00 

Hence it follows that u e VF 1 C (J) and further that (u, g) e D(A) and 

=i(«,g) = (f,g) . 
Further, we conclude by Fubini's theorem that 

r x 

dx 



c = -Nr 1 f 9(y)dy + HpUr 1 f P(x) \ \ X f(y)dy 

Ja Ja \_Ja 

= \\P\\i 1 \ p(x)dx f(y)dy- g{y)dy \ . 

yJa \_Jy J Ja J 

This implies that 

rX ( rb T rb ~| rb 

q(x) = - f(y)dy + \\p\\^ 1 \ p(x)dx f(y)dy- g(y)dy 

Ja {.Ja \_Jy J Ja 

for every x e I. Further, again by Fubini's theorem, it follows that 

rx rx rx 

u(x)=\ g(y)dy+\ p(y) q (y) dy + c p(y) dy 

Ja Ja Ja 
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rx rx r ru px 

g(y) dy + p{u) f(y) dy du + c p(y) dy 

Ja Ja \_Ja J Ja 

rx rx r rx 

g(y) dy + p(u) 

Ja Ja \_Jy 



f(y)dy 

rx ( rb r rb ~| rb 

+ \\P\\T l P{y)dy < p{x)dx f(y)dy- g(y)dy 

Ja (Ja \_Jy J Ja 

for every x s I. In addition, by Hoelder's inequality, we conclude that 

\u(x)\ < 2 (b - a?* [ \\Ph \\fh + Nb] < 2(6 - a)V2 ( i + \\p h) 
| g (x)|^(6-a) 1 /2[2||/|| 2 + ||p|| r 1 || 5 || 2 ] 

^2 (6 -a) 1 / 2 11^11^(1 + bli) ||(/,5)||2 
for every x e I. The last implies 

H| 2 <2(6-a)(l + |p|| 1 )||(/, 5 )|| , ||g|| 2 ^2(6-a)|b-||r 1 (l + |b1i)||(/,5)||2 
and 

\\(u,q)\\^2(b-a) MlHl + Mi) 2 Kf,9)h- 
As consequence, by ((L^.(/)) 2 — > (L^(/)) 2 , (/, g) — » g)), there is defined a 
compact bounded linear operator B. Since 

AB(f,g) = (f,g) 

for every (/, g) e (L 2 ^/)) 2 , the bijectivity of A follows as well as that A^ 1 = B. 
Finally, since Ar 1 is compact, A has a purely discrete spectrum. □ 

Corollary 6.0.24. Let a, b,p as in Theorem 6.0.23. Then U r (0), where 

r:=2- 1 (6-a)- 1 (1 + Ml) -2 . 
is contained in the resolvent set of A 
Proo/ For this, let A 6 U r (0). Then 

A- A = (1 - Ai _1 )i . 
By help of the previous Theorem 6.0.23, it follows that 

|A|-||i- 1 K|A|/r<l 
and hence that A — A is bijective. □ 
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Theorem 6.0.25. Let a, b e R such that a < b and f2 := J := (a, b). Further, let 
Pi,p2 e L°°(f})\{0} be a.e. positive and A±,A2 be the corresponding operators. 
Then 



Hr'-Vl^i 



IMIl 

Proof. Proceeds by direct calculation. 



||Pl||i + ||p: 



1 



2 1 



32 1 



\P2 ~ Pi ||l 



□ 



Corollary 6.0.26. Let a, 6 e R such that a < 6 and Q := / := (a, 6). Further, let 
Poo 6 L°°(f2)\{0} be a.e. positive. Let pi,p~2> ■ ■ ■ be a sequence of a.e. positive 
elements of L°°(n)\{0} such that 

lim \\p u - pa, ||i = . 

>oo 

In addition, let A±, A2, ■ ■ • be the associated sequence of self-adjoint linear opera- 
tors and be the self-adjoint linear operator associated to p^. Then 



lim || .A 

v— >co 



-1 



.4 



~ii 



. 



Proof. The statement is a simple consequence of Theorem 6.0.23. 



□ 



Corollary 6.0.27. Let a, b e R such that a < b and := I := (a,b) and / e 
Lq(I). Further, let poo 6 L°°(fi)\{0} be a.e. positive and pi,p~2, • • • be a sequence 
in £ such that 

lim \\p u - poo ||i = . 
v— >oo 

In addition, let Ai , A3 , . . . be the sequence of self-adjoint linear operators that 
is associated to px,p2, ■ ■ ■ and p v := l/p u for v e N*. Then A^^A^ , ■ ■ ■ 
and —p 1 DfA^ 1 , — p 2 Z?|A^ 1 , ... are convergent in L(L^(I), L^(I)) to 6,C e 
L(L^(I), L^(I)), respectively. In particular, B and C are given by 



(u) du 



f(y)dy 

■b r rb 



rx rx 

(Bf)(x) = p 

Ja \_Jy 

rx rO r ro 

+ IIpooIIi Pcc(y)dy poo(x)dx 

Ja Ja \_Jy 

(Cf)(x) = - [ f{y)dy + UPoollr 1 f ffpoo(x)dx 

Ja Ja \_Jy 



f(y)dy , 
f(y) dy 



for all x G / and every / 6 L^(I). 
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Proof. The statement is a simple consequence of Theorem 5.0.20 and Theorem 
6.0.23. □ 



7 Concluding remarks 

It is unclear whether results similar to those of the previous section can be ex- 
pected to hold in dimensions greater than 1. According to Theorem 5.0.17 and 
the subsequent example, and differently to the situation in one dimension, A is not 
injective when p vanishes on non-empty open subsets of the material. Hence there 
does not seem to be an obvious candidate for a limit of a sequence of A' 1 that is 
associated to a sequence in C approaching such p. Therefore, it is conceivable that 
such limits show a wider variety of phenomena than those in one dimension. This 
problem deserves further study. 

A final remark concerns the fact that it cannot be expected that general 'elliptic 
regularity theorems' hold for operators A corresponding to discontinuous diffusiv- 
ities p as a consequence of the condition that every element u from the domain of 
such operator satisfies pV w u 6 /J(Vq). For instance, the source function / in 
Example 4.0.1 1 is in W£(I) for every k e N, but u = A -1 f £ W£(I), where / is 
the open interval (—1,1) of R. 

8 Appendix 

In the following, proofs of the Lemmata 3.0.4, 3.0.6 from Section 3 are given. 
Lemma 8.0.28. (Partial integration) 

<f\d ek g\ = -{d ek f\g) 2 

for all (/, g) 6 Wq c (fi) x W^(Q) and k s N*, where e& denotes the k-th canonical 
unit vector of W 1 . 

Proof. For this, let k s N*. We define the sesquilinear form s : Wqq(Q) x 
W£(fi) -> Cby 

s(f,g):=(f\d e *g) 2 + (d e *f\g) 2 

for all (f,g) e Wq C (Q) x 1^(0). By the continuity of d ek , it follows the con- 
tinuity of s and by partial integration that s(f, g) = for all / 6 Qj°(f2, C) and 
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/ e C°°(fi,C) n W£(fi). Since C$>(Q,C) * ((C°°(fi,C) n is dense 

in Wo C (n) x W(j;(f2), this implies the vanishing of s and hence the validity of 
(3.0.4) for all (/, g) e W^Q) x W&Sl). □ 

Lemma 8.0.29. (Adjoints of gradient operators) 



(V ) = ^w\ w i c{n) , ( v «lwi c (n) J " v °* • 
Proof. Since Vo* is densely-defined, it follows that 

(v *r = Vo . 

For / 6 D( Vo ), there exists a sequence fi, f2, ■ ■ ■ in Cg°(f2, C) such that 
lim \\f u - f\\ 2 = , lim ||Vo/v - Vofhn = lim || VJ, - Vo/|| 2 ,n = . 

V—XX) 1^— >00 j;— >co 

Hence it follows that / 6 W 1 c (O) and that Vo/ = V w f. As a consequence, 

Vo c V^l^i^ • 

Further, for / e Wq C (Q), there is a sequence /x, /2, • • • in Cq 3 (Q, C) such that 
lim \\f v - f\\ 2 = , lim ||V /, - V w f\\ 2 = . 

!/— >CO I/— >00 

Hence it follows that 

(/,V„/)eG(Vb) . 

As a consequence, 

Finally, it follows the validity of (3.0.5, 1). The validity of (3.0.5, 2) is a simple 
consequence of (3.0.5, 1) and the closedness of Vo*. □ 
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